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1 Introduction 



1.1 General remarks 

The brightest representatives of art, say, Mozart, Mendelsohn, etc., even though 
they were known as wundcrkinds from early age, achieved their full maturity 
only at the time of their death, e.g. recall Requiem (for Mozart) or 9-th sym- 
phony for Beethoven, etc. In the case of Landau, his scientific abilities had 
become apparent very early and his international acclaim came to him when 
he was still in his 30ies. Nevertheless, like with other greatest, the most pro- 
found of his papers were written in a short period of only 3 years, from 1956 
to 1959 (e.g. see the complete list of Landau papers in this collection) prior to 
his tragic fatal car accident, which happened on 7th of January of 1962. As in 
the case of Mozart or Beethoven, these works stand in their own classes since, 
as if anticipating his own end, with these works he made an attempt to take 
a look into distant future trying to foresee contours of the ultimate " theory of 
everything". In his last paper "About Fundamental Problems", written in 1959 
on the occasion of Wolfgang Pauli's death, he summarized his own efforts as 
well as those whom he considered as the best from his generation of physicists. 
In retrospect, very much like Einstein who wrote his own obituary at the age of 
68, this paper happened to become Landau's obituary as well. 

About 50 years after this paper [1] was written, it is impossible not to be as- 
tonished by Landau's vision of the future, especially, in view of his own cautious 
words: 

"As is well known, in the case of strong interactions, theoretical physics at 
this moment is to a large extent powerless. Therefore, any attempts at making 
statements regarding future directions of this field should be considered as highly 
speculative (and, hence, prophetic) causing their authors to be in ever increasing 
danger of shooting straight into the sky." 

In spite of these cautious remarks, he goes on in this paper and makes 
"profecies" which happen to be higly accurate. What makes this paper truly 
remarkable is the fact that in it (in retrospect) he actually was speaking on 
behalf of the whole generation of physicists of his time. The famous photograph 
of Dau sitting next to Geo Gamow and Hcndrik Kramers at the Niels Bohr 
Institute seminar carries much more than the flavor of these exciting times. 
Looking back, one notices that the people sitting in the front row on this picture: 
Bohr, Heisenberg, Pauli, Gamov, Landau and Kramers, although of different age 
at the time when the picture was taken-all died within the same time span of 
about 10-15 years. Landau's paper is an obituary for W. Pauli (died in 1958), 
but then, N.Bohr died in 1962, H.Kramers-in 1952, G. Gamov -in 1966, Landau 
-in 1968 and only W. Heisenberg died later in 1976. These circumstances make 
the content of his last paper especially significant. This is so, because 50 years 
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later upon its reading, correctness of Landau prophecies is truly remarkable. In 
the rest of this paper we would like to explain why this is indeed so and, based 
on this, we also will make a few predictions of what lies ahead. 

1.2 Overview of the last Landau paper 

In his obituary, Landau was concerned about the future of high energy physics 
leaving aside predictions for other branches of physics. This is justifiable since 
historically, especially in Landau's time, the high energy physics was the ulti- 
mate source of inspiration for the rest of physics. Landau begins with discussion 
of utility of methods of renormalization and renormalization group. He notices 
that in the case of strong interactions use of the renormalization group meth- 
ods leads to paradoxical results: even if one begins with very large couplings 
at low energies, at high energies the coupling tends effectively to zero. This 
phenomenon of "asymptotic freedom" was discovered by Pomeranchuk in 1955. 
At the time of this discovery, to cure such a nulification the nonlocality of inter- 
actions was cautiously suggested. This nonlocality idea was not welcomed by 
Landau since, in his opinion, all consequences of the mathematical apparatus 
of quantum field theory made without any reliance on a specific model Hamil- 
tonian require locality and are supported by all available experimental data, in 
particular, by sucsessful use of dispersion relations. Furthermore, in support of 
his claim, Landau used the theory of random fragmentation developed by Fermi 
[2] for calculation of angular distribution of pions produced in (high energy) 
collisions of two nucleons. It required the application of concepts of statisti- 
cal thermodynamics to space-time volumes much smaller than those needed for 
development of the nonlocal theory. These ideas were checked experimentally 
with large degree of sucsesfH. Next, Landau notices that Heisenberg [3] recently 
expressed the opinion that the existing theory most likely should undergo sub- 
stantial changes caused by his S-matrix theory. These changes however should 
not be done at the expence of removing the locality requirement in view of the 
assumed validity of dispersion relations implying Lorentzian causality. Landau 
goes on and makes a cautious remark about uscfulnes of Hciscnbcrg's nonlin- 
ear field theory also discussed in [3j§. His reservations are based on results of 
his work with R.Peierls (done in 1930) which questions the existence of wave 
functions in relativistic quantum field theories. It should be said, however, that 
such negative result he attributes only to the theory of strong interactions. In- 
deed, he writes: " Operators ip, containing nonobservable information, should 
disappear from the theory and, since the Hamiltonian can be made only from 
ip s, we necessarily come to the conclusion that the Hamiltonian method for 
strong interactions outlived itself and should be hurried, of course, with all due 
respect". He goes on by saying that the " basis of new theory should be made out 
of diagrammatic technique, which uses only diagrams with "free" ends, that is 

2 Deatails will be discussed below. 

3 It should be nothed though that Landau avoids any explicit mentioning about this non- 
linear field theory and only remarks that Heisenberg's assumptions in his opinion are ques- 
tionable. 
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typical for scattering amplitudes and their analytical continuations. The physi- 
cal foundation of such a formalism is made of principles of unitarity and locality 
of interactions revealing themselves in analytical properties of fundamental ob- 
jects of this new theory, for instance, in all kinds of dispersion relations 

As result of such an approach, the old problem of " elementarity" of elementary 
particles looses its meaning since it cannot be formulated without interaction 
between particles." This point of view is in accord with that made by Heisen- 
berg in [3] . The above quotations taken from Landau paper comprise to a large 
extent its content. 

2 Impact of Landau paper on high energy physics 

Since the vaue of science for society is based on its ability to explain the present 
and to predict the future, we would like to apply these criteria to Landau's last 
paper. We begin with 

The theory of random fragmentation [2] . Development of this theory resulted 
in development of QCD and explanation of the phenomenon of asymptotic frce- 
We restore some details of this development following the recent review 
paper by Tannenbaum [4]. In view of this, for the sake of space, we squeeze the 
number of key references to the absolute minimum. In the early 60ies (the time 
when Landau was incapacitated but still alive) construction of proton acceler- 
ators with energies well above the threshold for anti-proton production, a zoo 
of new particles was discovered. Gell-Mann and Ne'eman noticed that particles 
sharing the same quantum numbers (spin, parity) follow the symmetry of the 
Lie group SU(3). This group is built of three generators which they associated 
with three fictitious particles (quarks): u- for "up", d-for "down" and s-for 
"strange"- all having spin 1/2 and fractional charge. According to the emerging 
picture mesons are made of two quarks while the barions-of three. This lead to 
the discovery of a new f2 particle in 1964 based on predictions of Gell-Mann- 
Ne'eman theory. This particle had an apparent problem with the Pauli principle 
since it was made of three s-quarks. To avoid this problem, it was suggested in 
the same year that quarks come in 3 colours. The major breakthrough came 
after it was realized that the SU(3) symmetry is NOT associated with three 
quarks but with three colours and that colour-charged gluons are the quanta of 
asymptotically free strong interactions which bind hadrons. 

QCD and the asymptotic freedom. Although the path to asymptotic free- 
dom begins with Pomcranchuk's paper written in 195^3, it had many twists 
and turns before it reached its final destination culminating in works by Gross 
and Wilczek [5] and Politzer [6]. The above mentioned breaktruth lead to re- 
newed attention to the Yang-Mills (Y-M) theory [7] which was untill that time 
in the shadows because of problems with its renormalizability. The major ad- 
vancements in solving the renormalization riddle for the Y-M fields occured in 

4 Much more about this theory will be said in Section 5. 
5 E.g. see Ref. [2] in Landau's paper [1]. 
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1967 (when Dau was still alive!). In this year Faddeev and Popov [8], follow- 
ing earlier general ideas by Dirac [9], came up with the correct perturbative 
scheme for the Y-M fields which takes into account gauge constraints and 
gauge redundancy. The famous Faddeev-Popov auxiliary fermionic fields (the 
Faddeev-Popov ghosts) were incorporated into gauge-invariant renormalization 
scheme independently by Slavnov in 1971 [10] and 't Hooft and Veltman in 1972 
[11]. In the mean-time using renormalization group methods Bjorken came up 
with his famous parton model in 1969 [12] which was sucsessfully tested experi- 
mentally and was instrumental in proving the asymptotic freedom for the Y-M 
fields [5,6]. In about the same period Veneziano, following Regge ideas about 
the likely form of asymptotics of scattering processes at very high energies, and 
in complete accord with Landau's predictions (regarding the central role of am- 
plitudes) came up with his famous amplitude in 1968 [13], the year of Landau's 
death. This amplitude, as is well known, gave birth to dual resonance models 
[14] and string theory [15]. Incidentally, much later on the whole diagrammatic 
machinery of perturbative calculations in quantum field theories (of whatever 
kind, in general, and of QCD in particular) was redone in the string-theoretic 
fashion [16] making practical calculations considerably easier. We shall discuss 
this subject in some detail in Section 4. In the meantime, we would like to 
explain why it was nesessary to come up with string theory in the first place. 

The standard model. The story begins with two names: Goldstone and 
Higgs. In 1961 Goldstone published his paper [17] on global symmetry breaking 
which is always being accompanied by the emergence of a massless particle, 
the "Golstone boson", while Higgs in 1964 published a paper on local symme- 
try breaking causing emergence of the massive particle subscqucstly called the 
Higgs boson [18]. In both cases, it is rather easy to recognize Landau's in- 
put since the global symmetry breaking phenomenon follows directly from his 
work on phcnomenological theory of phase transitions [19], while the breaking 
of local symmetry causing emergence of the massive Higgs boson is contained 
in his work (with V.Ginzburg) on the phenomenological theory of supercon- 
ductivity [20]. Breaking of local symmetry in this case allows to explain the 
Meissner effect in superconductors: intitially gauge-invariant Lagrangian of the 
Ginzburg-Landau theory of superconductivity (technically known as Lagrangian 
for scalar electrodynamics) formally looses its manifest gauge-invariance since 
the massless electromagnetic field acquires mass as a result of spontaneous local 
gauge symmetry breaking. Emergence of this mass causes the magnetic field 
to be expelled from the bulk of the superconductor. This is the microscopic 
cause of the macroscopically observed Meissner effect. Since many contributions 
to this centenary volume discuss this and (related to it) topic(s) in sufficient 
detail, we continue with other issues in our review. 

In their Nobel Prize winning lectures both Weinberg and Salam emphasized 
the importance of Goldstone and Higgs (and, consequently, Landau's) ideas in 
formulation of the unified theory of weak and electromagnetic interactions. In 
the Lagrangian of the Weinberg-Salam (W-S) model one can easily recognize 
the Ginzburg-Landau (superconducting) part responsible for the emerging Higgs 
boson. To make any serious calculations with this model requires use of the 
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rcnormalization group methods applicable to the non-Abelian Y-M theories in 
which gauge invariance is spontaneously broken. These methods were developed 
by 't Hooft and Veltman [11]. With help of these methods, it has become possi- 
ble to extend the W-S model in order to include strong interactions. The model 
which unifies weak, electromagnetic and strong interactions has become known 
as the standard model. Many predictions of this model were thoroughly tested 
experimentally. An excellent summary of both theoretical and experimental 
results related to the standard model can be found in a recent monograph by 
Bardin and Passarino [21]. More relaxed and up to date exposition of the 
satandard model can be found in [22]. 

Higgs boson and unified field theories. With all the sucsesses of the stan- 
dard model, it suffers from two major drawbacks: a) it needs Higgs (Ginzburg- 
Landau) field to be renormalizable, b) it does not account for gravity. Inclusion 
of gravity into existing unification scheme is complicated by the well known 
problems of renormalizability of gravity. Although, as was demonstrated by 
Utiyama [23], both gravity and the Y-M theory are obtainable as result of im- 
position of the requirement of local gauge invariance with respect to: a) external 
(Lorentz) group-for gravity or, b) with respect to internal (say, isotopic) sym- 
metries in the case of Y-M, on the underlying field-theoretic Lagrangian, from 
the point of view of perturbative renormalization group analysis, gravity escapes 
treatmets which are successful for the Y-M fields. As result, the string theory 
was proposed. It is a nonlocal theory, an option which Landau had anticipated 
but did not like, and is still in the process of development which, at the time 
of this writing, has lasted already for 40 years. More on this will be said in 
Sections 4 and 5. 

As we've mentioned already, the existing string-theoretic formalism can be 
adapted to the standard model [16] so that if the string theory is sucsessful, 
then the standard model is easily recoverable from it. Having said this, we 
silently assumed that the sucsess of string theory depends very much on its 
ability to incorporate and to treat gravity on equal footing with the rest of the 
forces. This leaves the door totally open for explanation of the existence of the 
Higgs boson. In the standard model, even though it is not yet discovered, its 
existence is not questioned^] because of the following: a) many predictions of the 
standard model not directly involving this boson happen to be correct [21], b) 
this boson is needed for this model to be renormalizable, c) only one such boson 
is sufficient for the whole machinery to work effectivelj0. The artificiality of the 
Higgs boson can be completely removed though if the already published results 
contained in our latest works [24,25] become commonly accepted. In these 
papers we demonstrated that the auxiliary Ginzburg-Landau (Higgs) functional 

6 Very much the same as the existence of gravitational waves, still to be discovered. It is 
interesting to notice that the Google search database has 599,000 entries on Higgs boson and 
513,000 entries on gravitational waves. In both cases one can find many claims that these 
objects will be discovered, if not today, then, at worst, by tomorrow morning. 

^Theoretically, nothing forbids us from studying the supersymmetric version of the Y-M 
model (and, hence, of QCD and the standard model). Such a model would require more than 
one Higgs boson for its renormalizability. Without discovery of the "standard" Higgs boson 
such a model thus far is only of mathematical interest (e.g. read the next section). 
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of the standard model can be (without any approximations !) rewritten as the 
Hilbcrt-Einsten functional for gravity. This means that such an equivalence 
converts the auxiliary functional for the Higgs field in the standard model 
into that for gravity thus making the already existing theory self-containecjf]. 
Such an approach differs substantially from that recently proposed by Lisi [26] 
since it does not require any additional efforts for developing the "theory of 
everything" . ft only changes our interpretation of the Higgs boson (very much 
like the work by Higgs which replaced more narrow treatment of the Ginzburg- 
Landau functional taylored for supercondactivity by mathematically broader 
tratment enabling to extract the underlying universality of the Higgs-Landau- 
Ginzburg phenomenon) . 

3 Impact of Landau paper on mathematical physics 

Exactly integrable low dimensional quantum systems. Although Landau was 
against the nonlinear field theory put forward by Heisenberg [3] and was advo- 
cating instead his diagrammatic methods for amplitudes, his intuition in this 
case was incorrect. Subsequent developments revealed that, actually, both (his 
and Heisenberg's) approaches are equivalent. Indeed, in 1979 the ground 
breaking paper by Zamolodchikov brothers was published in Annals of Physics 
[27]. In it, in complete accord with Heisenberg, the exact S-matrices were con- 
structed for a number of manifestly nonlinear relativistic models in 1+1 dimen- 
sions. Chronologically, this work can be considered as a culmination of efforts 
of many people, begining with McGuire [28], in 1964, then Yang [29], in 1967, 
then Baxter [30] in 1971, and many many others. The history of discovery of, 
now famous, Yang-Baxter (Y-B) equations is well documented in the collection 
of reprints assembled by Jimbo [31] where one can find as well many appli- 
cations of these equations to exactly integrable systems. In fact, the classical 
(quantum) system is exactly integrable only if the associated classical (quan- 
tum) Y-B equations can be solved [31,32]. The task of solving these equations 
was initiated by Belavin and Drinfel'd [33] and eulmunated in Drinfel'ds paper 
on Hopf algebras and quantum Y-B equation [34] . An excellet reference on Hopf 
algebras, quantum groups, Y-B equations, Drinfel'd associators and Knizhnik- 
Zamolodchikov (K-Z) equations is the book by Kassel [35] where one can also 
find connections of all these topics with theory of knots and links. The same 
topics are discussed in a somewhat more formal way (with emphasis on sym- 
plectic aspects of these issues not present in Kassel's book) in the monograph 
by Chari and Pressley [36]. A purely combinatorial, Heisenberg-style exposition 
of these topics can be found in our recent work on Heisenberg's honeycombs 
[37] in which our readers can find a very elementary introduction to K-Z equa- 
tions. These equations were discovered by Knizhnik and Zamolodchikov [38] 
in connection with their study of the Wess-Zumino-Novikov-Witten (WZNW) 
model (the nonlinear sigma model with the topological WZNW term [39]). Such 

8 Under such equivalence the Higgs boson should be renamed as the Higgs- Yamabe con- 
formon. 
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a model emerged as a by-product of developments in string theory and in Y- 
M theorj|^| but is of independent (from these theories) interest and value in 
view of its wide uses in condensed matter [39,42] and conformal field theories 
[43] where, incidentally, the Ginzburg-Landau description of such models exists 
[24,32,43]. In his lecture notes Varchenko argued [44] that all hypergeomet- 
ric functions of multiple arguments are solutions of the K-Z equations. Among 
these hypergeometric functions of multiple arguments are all Landau amplitudes 
as demonstrated by Golubeva [45] and will be discussed in Section 4. Hence, by 
studying these aplitudes as suggested by Landau, one inevitably runs into K-Z 
equations, quantum groups, Y-B equations, knots and links, etc. 

The Yang-Mills theory in various dimensions and Ginzburg-Landau model 
of superconductivity. The role of self-duality. Most likely, should he be alive, 
Landau might argue, based on the knowledge of his time, that the S-matrices 
obtained by Zamolodchikov brothers are of academic interest since they only 
involve nontrivial quantum field theories in one time and one space dimensions 
(which is reflected in his skepticism of the Heisenberg's program) . Such a con- 
clusion, however, is not correct because of the following. On one hand, these 
quantum models do have significance in real life in view of their experimentally 
supported connections with spin chains [39,42]. On another hand, attempts at 
quantization of higher dimensional nonlinear quantum fields are most likely to 
be in vain (in accord with Landau!). The combinatorial arguments in favour 
of this conclusion are presented in our recent works [37, 46, 47]. They also fol- 
low from Golubeva's paper [45], and will be further explained in the following 
section. In support of our point of view, we also would like to mention the 
problem of quantization of pure Y-M field. The nontriviality of such a quanti- 
zation already begins in 2 dimensions where serious study of pure Y-M living on 
Riemannian surfaces was initiated by Atiyah and Bott [48]. Their work caused 
an avalanche of subsequent works by mathematicians. In particular, in two 
dimensions, the complete solution of the pure Y-M field theory was found by 
Witten [49,50] whose work was significantly influenced by papers of Atiyah and 
Bott [48] and Migdal [51]. In dimensions higher than two the situation formally 
remains much less tractable since the two-dimcntional methods resist general- 
ization to higher dimensions. Nevertheless, in some special cases, e.g. when 
the Y-M is (anti)sclf-dual, the theory can be developed quite substantially. We 
believe, based on the results to be described immediately below, that only these 
realizations of the Y-M field make physical sense. 

We begin with very significant work by Belavin et al [52] done in 1975. In it, 
the instanton solution of the Euclidean version of the Y-M field theory was found 
using essentially ideas of self-duality. In the same paper the Chern-Simons (C-S) 
functional was used for the first time in physics literaturj^l. Shortly thereafter, 
't Hooft substantially improved calculations of Belavin et al (who essentially 
obtained only the leading saddle point result for pure Y-M) by making these cal- 

9 Where it originated from the Chern-Simons field theory [40,41] which itself is obtainable 
from the self-dual Y-M theory to be discussed below. 

10 The CS theory was used later on by Witten [40] in his Fields medal winning paper on the 
Jones polynomial in which many physical applications are discussed. 
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culations more realistic by incorporating fermions interacting via Y-M fields into 
these saddle point-type calculations [53].Atiyah, Hitchin, Drinfcld and Manin 
(AHDM instanton) developed results by Belavin et al much further using meth- 
ods of algebraic geometry. A very readable account of these and related works can 
be found in the review paper by Belavin [54]. Donaldson used essentially these 
AHDM results for development of his theory of topological classification of four 
manifolds [55]. In this theory the anti self-dual Y-M fields play prominent role 
since only such fields can be associated with complex structures on 4-manifolds. 
A very readable account of these developments including nice overview of the 
theory of 4 manifolds is given in a review by Soloviyov [56] . More up to date 
information, including that on Seiberg-Wittcn (S-W) theory, which enabled to 
simplify considerably the Donaldson theory is given in the paper by Tyrin [57]. 
The monograph by Nicolaescu provides an unprecedented amount of readable 
material relating the S-W and Donaldson theories [58]. At the same time, the 
S-W theory is just an elaboration on the phenomenological Ginzburg-Landau 
(G-L) theory of superconductivity as explained in detail in the monograph by 
Jost [590. In S-W theory the Abelian gauge field of G-L theory is replaced by 
the non- Abelian Y-M gauge field, the boson field is replaced by the spinor field 
and the G-L covariant derivative for scalar fields is replaced by that appropriate 
for the Dirac fields. 

In conclusion, we would like to mention that practically all known exactly 
integrable systems can be obtained from the solutions of the (anti) self-dual 
Y-M field equations of Donaldson's theory [62] . Similar results for the S-W 
theory are discussed in the monograph by Marshakov [63]. This means that 4- 
dimensional physics, including gravity (in this case, the self-dual gravity [64]) 
can be reduced to the study of two dimensional problems in the spirit of the 
program outlined by Heisenberg [3] and in accord with results of our latest works 
[37,46,470. 

Hitchin equations, Langlands duality, etc. With exception of the S-W theory 
which is an excellent tool for classification of 3 and 4 manifolds but, most likely, 
cannot be used in its original form for detailed comparison with experiment, we 
have not discussed the supersymmetric versions of Y-M theory. Mathematically, 
to study such versions is interesting and leads to issues such as mirror symmetry, 
Langlands duality, Hitchin's equations, etc. [65]. However, physically, such 
theories would require more than one Higgs boson. Since such a boson remains 
to be discovered even in the nonsupersymmetric casJ^l. there is no reason to 
go on in this review with description of ramifications caused by the effects of 
supersymmetry (may be, only with the exception of the Faddeev-Popov ghosts 
which do admit supersymmetric treatment [660 in the standard Y-M theory) . 

1 See also the books by Pismen [60], where emphasis on physical aspects of G-L theory is 
made, and by Yang [61], where mathematical aspects of [60] are discussed in detail. 
12 More on this is discussed in Sections 4 and 5. 

13 In spite of the very detailed experimental set ups and theoretical calculations backing 
these set ups. 

14 Known as the BRS-type symmetry named after Becchi, Rouet and Stora. 
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4 Impact of Landau paper on mathematics 

4.1 General remarks 



In his paper [1], Landau emphasized the most the results of his last long and 
deep paper on analytical properties of the vertex parts (essentially, the scattering 
amplitudes) in quantum field theory [67]. This work can also rightfully to be 
considered as Landau's last paper. In the review by Golubeva [45] our readers 
can find a nice summary of the impact of this paper on mathematics up to 1976. 
Accordingly, to avoid duplications, we shall discuss in this section developments 
which either happened after this date or which are not discussed in Golubeva's 
paper, but arc of importance. 

In particular, we would like to mention some inconsistency in Landau's rea- 
soning. On one hand, he wrote his paper [67] in order to make some progress 
in the theory of strong interactions. On another, the mathematical methods 
based on analysis of Feynman's diagrams which he was employing does not dis- 
tinguish between "strong" and "weak interactions. Because of this, we shall 
discuss analytical properties of Feynman's vertex diagrams for any quantum 
field theory. Using Golubeva's results we are forced to admit that, whatever 
this diagram might be, there is an equation for a hypergeometric-type function 
of multiple arguments whose solution is given by the vertex part in question. 
Varchenko demonstrated [44] that all solutions of hypergeometric equations of 
multiple arguments are solutions of the corresponding Knizhnik-Zamolodchikov 
equations. These solutions are all also expressable in the form the Aomoto- 
Gclfand type hypergeometric integrals [68-70] associated with GKZ equations 
(also of hypergeometric type) to be discussed below. Knowledge of analytical 
properties of the vetex part allows us to connect it via the optical theorem with 
the experimentally observable cross section [71]. Hence, cither such experimen- 
tal data can be used to reconstruct the vertex part via Kramers-Kronig-type 
dispersion relations, or the theoretically obtained vertex part (as some solution 
of the K-Z or GKZ equation) can be used for making predictions about the 
scattering crossection. In both cases only 2 dimensional results are actually 
used in accord with remarks made in the preceding section. The question then 
emerges: is there a way to classify all solutions of K-Z equations so that all 
meaningful scattering amplitudes can be obtained, recognized and classified? 
In part, such a classification might follow from works by Bclavin and Drinfcl'd 
[33] and Reshctikhin and Wiegmann [72]. And, because of this, all quantum me- 
chanical and field-theoretic models can be treated by the group-theoretic (which 
are in the essense topological [35]) methods involving theory of knots and liks. 
Such point of view is developed in our recent paper on Heisenberg honeycombs 
[37] and has found its implcmctnation in the field-theory case in the works by 
Connes and Kreimer [73,74], Connes and Marcolli [75], Kreimer and collabora- 
tors [76], etc. We would like to reobtain some of these results having in mind 
their applications in the next section. 

Before doing so, it is useful to remind our readers about developments in 
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mathematics, which were taking place exactly at the time when Landau's pa- 
per [67] was written. Because of these developments, his paper was immediately 
appreciated by mathematicians and made a big impact on mathematics. In 
1958, just less than a year before Landau's paper was written, Rene Thorn was 
awarded the highest award in mathematics-the Fields Medal-for his contribu- 
tions to the theory of cobordism [77] . This theory had played an important role 
in S-W and Donaldson theories, as well as those involving Jones polynomial [40], 
etc. Nevertheless, its is the theory of singularities and catastrophes developed 
by Thorn between 1954-1962 which is immediately relevant to Landau's work 
[67] . That this is indeed the case can be seen from the monograph by Pham [78] 
published in 1967. It is essentially the enlarged version of his earlier publication 
[79], which, in turn, is part of a still earlier paper by Fotiadi, Froisaart, Las- 
coux and Pham published in the top mathematical journal "Topology" [80]-all 
devoted to Landau singularities. Mathemaically, all these works are based on 
results by Thorn, Leray, Milnor and Lefshetz done in the time period around 
publication of Landau's paper. A systematic exposition of these mathematical 
results can be found, in addition to works already cited, in the monograph by 
Arnol'd, Varchenko and Gussein-Zade [81]. In this book, in addition, one can 
find the information about the so called mixed Hodge structures and period 
mapping-results of contributions of Dcligne and Griffiths. An excellent introduc- 
tion to these topics can be found in the monograph by Carlson, Muller-Stach 
and Peters [82]. All these works were employed in our analysis of analytical 
properties of Veneziano and Veneziano-Like amplitudes [83-85] eventually cul- 
minating in [47] replacing more traditional string-theoretic models [15] by the 
tachyon-free model existing in the usual space-time dimensions. We mention 
this in view of one of the latest papers by Kreimer and Bloch [76] in which 
they noticed relevance of the mixed Hodge structures to more traditional field 
theories requiring renormalization. Since Veneziano postulated (guessed) his 
amplitudes in 1968 [13] many attempts were made to derive these amplitudes 
field-theoretically by studying high energy asymptotics of scatterring processes 
where, as is commonly believed, the scattering amplitudes should exibit the 
Regge-type behaviour [86,87]. Accordingly, we would like to demonstrate now 
that, in fact, the Veneziano amplitudea 15 ! are the backbones of scattering ampli- 
tudes of every imaginable field theory. This fact will be used in the next section 
in which we are going to demonstrate much wider uses of Veneziano amplitudes 
than just those in the scattering processes of high energy physics. 

15 We shall call scattering amplitudes as "Veneziano amplitudes" since they all can be 
expressed in terms of, may be, linear combinations of Euler's beta functions of multiple ar- 
guments. The true Veneziano amplitude has in addition the Regge-type parametrization of 
these arguments. 
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4.2 Veneziano amplitudes and Feynman diagrams 

We begin with the key identity attributed to Feynman [88], page 83, 



1 



r(A ) 



n 



ii n n*i) ' 

3=1 3=1 



«(i-E«i)n«? j_1 



3=1 



3 = 1 



E 

v^ 1 



(4.1) 



Here Aq = E A? an< ^ ^X^j) ^ s Euler gamma function. The factors are given, 



as in Landau's paper [67], i.e. as Aj 



- . In the present case a slightly more 



general identity is used 16 !. Such an identity is employed for whatever Feynman 
integral of the type 

Bd d Pl d d p 2 ■ ■■ 



I 



with B, in the case of fermions, being some polynomial in k[s with ki being 
some momenta related to a given line on the diagram. Clearly, k% may be a 
combination of both the internal momenta, i.e. pi, and the external momenta, 
e.g. so that a given diagram is some function of the external momenta and 
dimensionality d of space (or space-time). Use of (4.1) in (4.2) significantly 
simplifies calculations. Details are discussed below. 

Looking at (4.1) we notice that in the case if Ay = 1 Vj we obtain an 
important identity 



r(A ) 



n 



II r(A,j 3=1 

3 = 1 



da, 



fa-E^II^" 1 - 1 



(4.3) 



3=1 



3=1 



This identity defines the Poisson-Dirichlct (P-D) probability measure to be dis- 
cussed in some detail in Section 5. Furthermore, without the normalization 



factor 



r(A ) 



-, we recognize in the obtained expression the Veneziano ampli- 



n r(A, 

3 = 1 

tude [83-85]. From this observation the following conclusion can be drawn: 
All Feynman's diagrams containing loops mathematically are averages of the 
stochastic processes of the Poisson-Dirichlet tvp(P^I.Because of this, all such di- 
agrams can be rewritten as linear combinations of Veneziano amplitudes. We 
would like to provide sufficient evidence that this is indeed the case. This task 
is greatly simplified by the significant progress made to date in actual calcu- 
lations of Feynman's diagrams. Recent papers by Bogner and Weinzierl [89] 
and Weizierl [90] contain a nice summary of these efforts. To these papers, we 



16 In Landau's paper all X^s are equal to one. 

17 In Section 5 we shall demonstrate that the same is true for traditional qauantum mechanics 
and many other disciplines. 
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would like to add the paper by Tarasov [91] on hypergeometric representation of 
a special class of Feynman's diagrams as well as papers by Smirnov and Smirnov 
[92] on the use of Grobner bases for calculation of Feynman integrals. Reading 
these papers and [16] provides a broad panorama of efforts made in practical 
calculations of Feynman diagrams to date, very much in accord with and in 
the spirit of Landau's predictions. 

For the sake of space, and without looss of generality we would like to con- 
sider only the bosonic-type diagrams (that is, the diagrams for which the factor 
B in the numerator of (4.2) is identically equal to one). This is totally justi- 
fied [89] since the diagrams which contain nontrivial _B-factors can be reduced 
to those in which B = l.In view of this, we are left with the calculation of 
integrals of the type 

I(a;s,t,u,...) = J (4.4) 

{P> A >) 

where a={ai, ...,«„} and s, t, u, ..are some kinematic invariants made of exter- 
nal momenta. Calculation of such integrals can be found in [89] and, therefore, 
are of no immediate concern to us. We shall consider their analytical properties 
irresppective to a particular outcome of computation of I(a;s,t,u, ...).To do so, 
we need to employ a few facts from the theory of functions of several com- 
plex variables. This is needed for justification of the Laurent-type expansion of 
I(a;s,t,u, ...).If such a justification is found, then we can represent I(a;s,t,u, ...) 
as 

oo 

I(a;s,t,u, ...) = ^ ci(s,t,u...)a l (4.5) 

/ — — oo 

with 1= [h, l 2 , .. .,/„]; a 1 = c^ 1 • •• a 1 ™, h + h + ■• + l n = h ci(s,t,u...) = 
c h,l2,...,l n ( s >ti u > ■■■)■ Clearly, in actual computations we do not anticipate that 
we shall use an infinity of coeffcients c/. But, whatever they are, if such ex- 
pansions do exist, then any Feynman integral is a Poisson-Dirichlct average 
< I(oc;s,t,u, ...) > with the probability density defined in (4.3). Hence, the 
task now lies in providing a justification that, indeed, such an expansion does 
exist. 

To do so, we need to recall a few facts from the theory of functions of many 
complex variables. In particular, we need to recall some models of complex pro- 
jective space CP™. These are discussed, for instance, in [93] and were employed 
in our works [83-85] on Veneziano amplitudes. 

Let C™ +1 be a complex n + 1 dimensional space and let u> — (loq, ...,ui n ) 
be a point in such a space. If we identify two arbitrary points to' and ui" in 
C n+1 via the relation w' = Aw", where A is some nonzero complex number, 
then the whole space C n+1 can be subdivided into equivalence classes with 
respect to previously mentioned relation. A quotient C n+1 /(ui' = Xui") is the 
complex projective space CP™. Each complex line in C™ +1 can be characterized 
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by the unit vector oj° = - — r . Every line made this way in C" +1 is a point in 

M 

CP™. Because of this, we notice that points in CP™ are also points in C™ +1 . 
This fact can be used for construction of another model representing CP™. 
For this purpose, let the line in C" +1 be given parametrically as z v = iv (, 
v = 0, ...,n with the parameter £ to be determined as follows. If an equation 

n 

for a complex sphere S n in C™ +1 is given by ^ z u z u — 1, then, we obtain: 

2 ™ 2 2 

ICI S \ w v\ = 1> implying ICI = 1 which is just an equation for a circle S. 
This means that CP™ can be determined as quotient of S n /S. Because of this, 

n 

we can define the deformation retract of CP™ as = 1 by identifying ti 

with |w? | .Thus obtained equation is an equation for a n— simplex A. We went 
into these details only because the 5 constaint in the P-D measure (4.3) tells 
us that the integration in (4.3) is done over the simplex A and, hence, such an 
integral lies, in fact, in the complex projective space (in the case of (4.3), in 
CP™ -1 ) in accord with Golubeva's paper [45]. Because of this, the meaning 
of auxiliary a variables in (4.5) changes from being real to becoming complex. 
This observation provides the needed justification for our use of the theory of 
functions of many complex variables. After recognition of this fact, we have 
to write down the correct differential form for CP™. To do so, we use results 
of Griffiths discussed both in Golubeva's paper [45] and in our paper [83] on 
Veneziano amplitudes, e.g. see Section 3.1. of [83]. Thus, we obtain, 

n 

w (w) = J2(-iy Wl dw[i], (4.6) 

i=0 

where dw[i] — dw\ A • • • hdwi-\ A dwi+i A • • • A dw n . Therefore, up to a constant, 
any nontrivial Feynman integral involving loops is expressable as 

- n 

I(s,t,u,..)= I(a;8,t,u,...)Y[ay- 1 «i(a) (4.7) 
J r j=i 

with r being some homology cycle in CP™ (very much like in the more familiar z- 
plane case one considers an integral of function f(z) along some contour which 
does not cross the singularities of f(z)) to be determined shortly. For this 
purpose, we notice that if the Laurent expansion (4.5) can be used in (4.7), then 
we end up with the sum of terms each of which is the Veneziano amplitude. 
From our work [83] we know than that each Veneziano amplitude represents 
some period on the Fermat hypersurface in CP" 

T(N) :a# + ...+a£ = 0, N = l,2,... (4.8) 

so that r represents one of the periods on such a surface (very much like in the 
case of the complex torus T there are two periods). This fact connects the theory 



14 



of such integrals with the period mapping theory, Gauss-Manin connections, 
hypergeometric functions of many variables, etc, e.g. see our paper [83] and the 
original works already mentioned, e.g. [68,69,81,82]. Such a conclusion depends 
on the existence of the Laurent-type expansions of the type given by (4.5) (used 
in (4.7)). Fortunately, the existence of such a result also follows, albeit implicitly, 
from the Lemma and Theorem in Section 1.4 of Golubeva's paper [45]. 

We would like now to discuss this topic in some detail. In particular, fol- 
lowing Shabat [93], we notice that upon rescaling, the differential form iv(w) 
behaves as follows 

*(«,) = ( -ir-i< + M-)A-Ad ( ^) Ad ( A- Ad ( ^) , (4.9) 

implying that a combination u(w)/F(w), with F(w) being some homogenous 
polynomial of degree n + 1, is scale- invariant. This is the result cited both in 
our work [83] and in the paper by Golubeva [45]. It belongs to Griffiths [94] 
who proved that in the projective space CP™ any closed differential form should 
look like $ = with degrees of homogenous polynomials P and F related to 
each other as degP + n + 1 = degF. Feynman diagrams discussed in Section 
1.6 of Golubeva's paper are all scale- invariant and obey Corollary 2.11. of 
Griffiths paper [94]. Scale invariance implies renormalizability and vice versa. In 
her paper Golubeva talks about convergence of Feynman integrals as a function 
of their dimensionality without paying attention to the renormalizability issues. 
Hence, contrary to Landau expectations, the diagrammatic technique he was 
advocating is allowing us to look only at the renormalizable theories. This is 
not a significant drawback though since, as we mentioned already, the Landau 
program also brings to light the hypergeometric functions of multiple arguments 
(to be futher discussed below) and, with them, the Knizhnik-Zamolodchikov 
equations [38, 44] and, hence, the WZNW model, etc. The scale-independent 
integrals are periods of some varieties, e.g. of the Fermat-type living in the 
complex projective space. In accord with general methods developed by Griffiths 
and nicely summarized in [82], it is of interest to study these periods as functions 
of parameters, in the present case-the parameters are the kinematic invariants. 
In addition, as in the theory of ordinary integrals calculated by the method of 
residues, in the multidimensional case there is analogous theory developed by 
Leray [79,80,93] and is also discussed in Golubeva's paper. But, as we know 
from the simpler, one-dimensional case, use of the methods of residue theory 
is essentially equivalent to our ability to obtain the Laurent expansion for a 
function in question. Hence, we come back to the expansion (4.5). 

In the case of the theory of functions of one complex variable the procedure 
for obtaining the Laurent expansion is well known. Such a procedure is not 
immediately transferable to the case of many complex variables though. To do 
this, one need to ask a question: what is the analog of the Cauchy formula in 
the case of many complex variables? Surprisingly, there are many analogs of 
this formula in the multivariable case. Let us recall what is involved in the 
one-dimensional case. For a domain D with boundary dD in z-plane and the 
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well behaved function f(z) in this domain we have the Cauchy formula 



dD 



It can be used only for z G D. For z outside D the result is zero since the function 
under the integral is holomorphic. With help of (4.10) the n— th derivative of 
f(z) is obtained as 



f(n)( z ) = 2!L I WW (4 11) 

dD 

This result is used for both the Taylor and Laurent expansions of f(z). By 
complementarity principle we expect that the same holds true in the multidi- 
mensional case. Hence, we have to find a multidimensional analog of (4.10) 
first. The requirements of scale invariance of integrals of the type given by (4.7) 
formally leave us with not too much choice. Our experience with Veneziano 
amplitudes [83] tells us, however, that the situation is considerably trickier. 
Indeed, taking into account the identity 

B(x, y)B(x+y, z)B(x+y+z, u)-B(x+y+-+t, I) = IgM— ll^ , (4.12) 

e.g. see equation (3.28) in [83], one can notice that study of analytical properties 
of any Veneziano amplitude can be reduced to that of the elementary Euler's 
beta function B(x,y) which is just a period of the simplest Fermat variety 
±Zq + Zi + z^ = so that the period integral in CP 2 representing this beta 
function can be written as 

z i lz 2 2z o° ( dz i dz 2 dz dz 2 dz dzi 
± z n + z n + z n( Zl a Z2 Zq a Z2 + ZQ A zi j, V-U) 

e.g. see (3.10b) of [83]. Here ci,C2 and Co are the same as x and y in B(x,y) 
while an auxiliary parameter cq is used only in the projective form of the beta 
integral. In actual comutations one should eventually make a transformation to 
the affinc form where one has to put zq = 1 so that the final result is independent 
of Co- Deatails can be found in [83]. For this (projective) period to make sence 
mathematically the following arguments should be applied. In the integral one 
should make the following replacements: Zk — ► z^ 3 , where f = exp(±i-^) 
(1 < 3 < N — 1) and k = 0, 1, 2. Substitution of such an ansatz into (4.13) and 
requiring that the integral / be independent of £ J leads to the (Veneziano-type) 
condition 

jici + J2C2 + joco = N (4-14) 

with j' k s being in the range specified above. Such a condition makes the period 
I nonsingular. Physically, however, this case is uninteresting since the nonsigu- 
lar expressions cannot be used for scattering amplitudes. The poles in such 
amplitudes (that is, the resonances) are being used in the optical theorem for 
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calculation of experimentally observed cross sections, e.g. read [71]. The way 
out of this apparent difficulty is explained in our work [83] , see also our related 
work [95] . In view of this, some caution should be exercised when one wants 
to use the Laurent expansion (4.5) in (4.7). This observation shifts attention 
from actual computations of coefficients in the Laurent expansion (4.5) towards 
studying of general aspects of such calculations. Our task in this section lies in 
proving that any scattering amplitude is a linear combination of the Veneziano- 
like amplitudes. This proof is nonconstructive (that is, it does not allow the 
effective computation of amplitudes) but it is essential for the discussion pre- 
sented in the next section. Again, we would like to remind our readers that what 
we call "Veneziano amplitudes" are in fact Euler's beta functions of multiple 
arguments. The true Veneziano amplitudes involve Regge-type parametrization 
of these arguments. Analysis made in [83] indicates that in the case when the 
Laurent expansion (4.5) has countable infinity of terms with negative powers, 
this counable infinity can be squeezed into just one (truly Veneziano) amplitude 
using the appropriate Regge-type parametrization. Such a situation is typical 
for the meson and hadron physics. In the case when such an expansion has just 
a few terms with negative powers, the situation should be typical for quantum 
mechanics (see next section), quantum electrodynamics and weak interactions. 

After these general remarks, finally, we are ready to provide needed (non- 
constructive) proof. For this purpose we need only to write the multidimen- 
soional analogs of the known Cauchy formulas (4.10), (4.11). In literature there 
are three types of multidimensional Cauchy-likc formulas-all reduceable to 
(4.10), e.g. see [93]. We prefer to use the Martinelli-Bochner (M-B) formul£0 
It is given by 

/(z) = J /(CWb(C-z) (4.15) 

dD 

with the M-B form ojmb being defined as 

"mb(C - z) = -j^-yr ^ ^— ^ (4.16) 



z 



and, as usual, \z\ 2 = ■ For n = 1, because of this, we obtain 



1 C-z 1 1_ 



lomb((-z) = —^ 7z^=——^, (4.17) 



as required. The analog of (4.11) in the present case was obtained by Andreotti 
and Norguet (e.g. see [93], page 229) and is given by 

/( k >(z) = j /(CH(C-z) (4-18) 

3D 



18 The more general Cauchy-Fantappie formula discovered by Leray in 1956 is reducible to 
that of M-B in the special case w = z, e.g. read page 158 (after equation 10) of [93]. We are 
interested just in this special case in this work. 
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with the form Wk being given by 

k! (n - 1)! ^ (-!)*& - ^) fc - +1 rfC" +I H A d( 

" k (C ~ z) = (2^f ^ ' (4J9) 

where k = {k%, k n } , a = {a±, ...,a n },I={l,...,l} and k! =k\\ ■ ■ ■ k n \. 

Thus, the Laurent expansion (4.5) does exists so that all scattering processes 
in quantum field and string theory are obtainable as averages of the stochastic 
Poisson-Dirichlet type processes. These are going to be discussed in some detail 
in the next section. Before doing so, we would like to conclude this section with 
a brief discussion of hypergeometric equations of multiple arguments associated 
with Feynman diagrams in order to bring Golubeva's paper [45] up to date. 

For the warm up, let us consider a calculation of the following generic integral 

j w = ^Jm d(: (4 - 20) 

7 

where P(z) is some poplynomial: P(z) = z n + a 1 z™~ 1 H ha„ and the contour 

7 is chosen in such a way that all zeros of P(z) lie inside. Evidently, 

P(±)I( Z )=Q. (4.21) 

az 

This equation should be supplemented withthe boudary conditions. They come 
from consideration of integrals of the type 

'"Mi— S?/j£j* ("») 

7 

Since k < n, the function pT^s nas a ^ the point z — oo zero of order n — k. 
For k < n — 2 we then obtain res^-oo — 7-7 = 0. While for k = n — 1 we obtain, 

c_ P(C) 

^-Tl— 1 ^-Tl— 1 

-— y ~ I for C ->■ 00. Hence, res c=oc -^y = -1 Therefore, /("-^(z) | 2=0 = 1 

and the rest of derivatives are being zero. This generic example can be broadly 
generalized. This task is accomplished in a series of papers by Gelfand, Kapra- 
nov and Zelevinsky (GKZ) whose results are summarized in [96]. In addition 
to this reference, we shall follow in part more relaxed exposition [97] describing 
results of GK2I 19 I. For the sake of space, we are not discussing the related to 
GKZ issues of mirror symmetry. A very accessible exposition of this topic is 
given in our work [95] to which we refer our readers for details. Reading of 
this reference should be sufficient for understanding of finer details of mirror 
symmetry presented in [97, 98]. 



19 In our exposition we fill in many gaps in presentation which one might encounter while 
trying to read [97] . 



18 



Although not mentioned in the review [96] , the original motivation for sudy- 
ing the GKZ hypergeometric functions came from observation made by GKZ in 
the earlier paper [98] that the integrals of the type 

F a {a,f3;P)= f ]J P;(z 1; x k )^x{ 1 • ■ ■ x{ k dx ± ■ ■ ■ dx k (4.23) 

i 

with a = («!,... ,a m ) e C m ,(3 = (ft, f3 k ) G C fc , P(x) = where 
x w = a;" 1 ■ • ■ x£ fc so that w = (wi, Wfc) 6 Z fc and a is some carefully chosen 
k-cycle (whose construction is described in [96, 98])-all are integrals of the quan- 
tum field theory (qft) already discussed in Landau's papers [1, 670 Although 
GKZ promised to write a separate paper devoted to the qft integrals only, to 
our knowledge such paper was never written. Evidently, our equation (4.7) falls 
into the category of GKZ integrals as required so not much else can be said. 
To deal with integrals of the type given in (4.23) consider the following integral 
very similar to our (4.7). It is given by 

4 m) (u) = / P u (x) m ^i • • • ^ (4.24) 
J (j Xl x k 

Here m E Z,P„(x) = M aX a , a=(ai, ...,& k ),A = {ai,...,ajy} is some finite 

subset of Z fe (usually being comprized of vertices of some convex polytope V 
thus leading eventually to the mirror symmetry arguments as explained in [97, 
98]). For simplicity a = &\ x • ■ -a k is a product of k circles Ui,i = 1 — k, 
<Ji E C, each centered at so that P u (x) ^ V (xi, ...,£&) S a\ X • ■ -a k . By 
differentiation under the integral sign we obtain 

^V) _ „ f x ap u(xr -i^l . . . (4.25) 



dui J a xi x k 

To develop this result further some knowledge of solid state physics is helpful. 
We would like to remind our readers of some useful facts to help the under- 
standing of what follows. In the case of more familiar 3 dimensional lattices 
one can choose (rather arbitrarily) a cell made of vectors ei,e2 and e3 so that 
any vector A of such a 3d lattice is decomposable as A=niei+n2e2+n3e3 with 
li; G Z. Suppose now that we translate this lattice as a whole by vector b. 
This will cause us to write: A=niei+n2e2+n3e3 + b(ni+n2+na). To make 
such a vector representation well defined we mayrequire that ni+n2+n3 = 0. 
If, in addition we fix the origin, this would result in an additional constraint 
(determining the location of the origin) niei+n2e2+n3e3 = 0. In the present 
case we have a lattice made of basis set {ai, ...,ajy} = A, so that the above two 
conditions are translated into 

d H h l N = (4.26a) 



20 In mathematics these integrals have become known as Aomoto-Gelfand-type integrals as 
mentioned already. 
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and 



[iai H h l„a N = 0. 



(4.26b) 



Since the vector a in (4.25) belongs to the set A, the condition (4.26b), when 
applied to the integral (4.25), leads to the following GKZ equation 

This equation is nesessary but is not sufficient for description of the equation 
for hypergeometric function of multiple arguments since it does not take into 
account the homogeneity of such a function. To account for homogeneity we 
need to consider what will happen to the integral li m \u) upon rescaling of 
the external parameters u. For this purpose some results from our papers [84, 
85] summarized in [100] are helpful. Omitting technicalities, we would like to 
present here only the "bootom line". To do so, we would like to consider again 
the Laurent polynomials of the type 

/(x) - £ Ua x a . (4.28) 

Here S a is the polyhedral cone associated with the Newton's polytope V . In 
connection with thus defined /(x), consider now a general definition of a quasi- 
homogenous function of degree d with exponents j 1 , ...7„. Such a function is 
defined by 

f(s^x 1 ,...,s^x n ) = s d f(x 1 ,...,x n ). (4.29) 

If we apply such a requirement to the individual term of the Laurent polynomial 
in (4.28), we obtain 

n 

X/',",, J = 1 - N - ( 4 - 3 °) 

where N is the number of terms in the Laurent expansion (4.28). Equation (4.30) 
is an equation for the hyperplane. Different hyperplanes may have different d'jS. 
The polytope V is made of a collection of N hyperplanes. Finally, let us 
differentiate (4.29) over s and let s = 1 in the end. Thus, we obtain 

i 

At this point our readers already can correctly guess that the required supple- 
mental equation(s) to (4.27) are equations describing the polytope. In view of 
(4.30), (4.31), these are k equations of the type 

[fa j tt j i-_d]4")(„)=0. (4.32) 
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We would like to illustrate these general facts by more familiar example of a 
usual hypergeometric function. For this, we need to choose the A-system of 
basis vectors as follows 



1 




-1 












1 









1 







1 














1 



(4.33) 



while for the d- vector we choose d=(l-c,-a,-b). Under such circumstances, the 
GKZ equations can be explicitly written as 



d 2 



d 2 



du\du2 du^dui 
d d 



$ = 



mi o u 2 ^ — 

OUi OU2 



Ui 



d 


d 1 


dui 


ou 3 _ 


d 


d 1 




ou 4 



$ = (1 - c)*, 
$ = -a$, 



From the second equation we obtain 

d 2 



d 2 d 
., $ = u^iui— - +c— )$ 
OU1OU2 au\ oui 



From the third and fourth equations we obtain as well 

d 2 



du^dui 

Finally, in view of (4.34) we obtain 



d_ 

dm 



d_ 

dm 



(4.34) 

(4.35) 
(4.36) 
(4.37) 

(4.38) 
(4.39) 



U 3~ 1 «I 1 ("l^2 +(! + « + b ) U ^ + ab )^ = ^^Qu 2 + C d^^- ( 4 ' 4 °) 

In this equation we have to set = U3 = 114 = 1 and u\ — z in order to get the 
familiar equation for the Gauss hypergeometric function 



d 2 



z(z - 1)^ + [(a + b+ l)z -c] + ab 



f(a,b,c;z) = 0. (4.41) 
Recall now that the integral representation for such a function is given by [101] 



f(a,b, c;z) 



T(c) 



r(&)r(c-&) 



c-b-l 



tz) a 



=< (1 -tz)- 1 >, (4.42) 
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where the Poisson-Diriclet average < • • • > is defined by 



1 



<•••> = 



r(c) 



/ 



c-b-l 



(4.43) 



T(b)T(c-b) 



o 



in accord with (4.3). By the principle of complementarity this means that all 
Feynman's diagrams, including those for vertex parts determining the scattering 
amplitudes are the P-D averages involving finite (e.g for a < in (4.42))or 
infinite (e.g. for a > in (4.42)) combinations of Veneziano amplitudes in the 
sence already described. Other implications of the results we have just obtained 
are discussed in the next section. 

5 Outlook: Impact of Landau's last paper on 
sciences in new millenium 

5.1 General comments 

Although Veneziano had guessed (postulated) his amplitude [13], and some au- 
thors have criticized such an approach to scattering processes of high energy 
physics [102], the results of the previous section demonstrate that Veneziano 
amplitudes or their linear combinations are intrinsic objects of high energy scat- 
tering processes. The question arises: If this is the case, what else can be said 
about these amplitudes? We noticed already in (4.3) that all such amplitudes 
are P-D averages. Hence, the task now lies in discussing generalities of such 
types of stochastic processes. Very fortunately, without any reference to high en- 
ergy physics this task was to a large extent accomplished. References [103-105] 
provide an excellent introduction into the theory of the P-D processes which 
play the central role in the theory of random fragmentation and coagulation 
processes. It is sufficient to type " Poisson-Dirichlet" using Google search engine 
in order to find about 52,000 entries. Such an abundance of entries is caused 
by the fact that many disciplines from the theory of spin glasses to computer 
science, from linguistic to forensic science, from economics to population genet- 
ics, from chemical kinetics to random matrix theory, etc. -all involve the P-D 
distributions. Remarkably, none of these references mention applications to the 
high energy physics, quantum mechanics or qft. The role of the coagulation 
fragmentation processes in high energy physics was recognized some time ago 
by Mekjian, e.g. see [106] and references therein. His works do not contain, 
however, arguments and results presented above, in Section 4, and therefore can 
be considered as complementary to ours. We shall say more about this below, 
in this section. 



22 



5.2 Random fragmentation and coagulation processes and 
the Dirichlet distribution 



We begin with some known facts from probability theory. For instance, we 
recall that the stationary Maxwell distribution for velocities of particles in a 
gas is just of Gaussian-type. It can be obtained as the stationary solution of 
Boltzmann's dynamical equation maximizing Boltzmann's-type entropjl^]. The 
question arises: Is it possible to find (discrete or continuous) dynamical equa- 
tions which will provide known probability laws as stable stationary solutions? 
This task will involve finding of dynamical equations along with the corre- 
sponding Boltzmann-like entropies which will reach their maxima at respective 
equilibria for these dynamical equations. We are certainly not in the position 
in this closing section of our paper to discuss this problem in full generality. 
Instead, we focus our attention only on processes wihich are described by the so 
called Dirichlet distributions. These originate from the integral (equation (2.8) 
in our work [83] on Veneziano amplitudes) attributed to Dirichlet, that is 

V(x u ...,X n+1 ) = [■■•[ < 1_1 ---< n-1 (l-Ul Un^+^dUf-dUn. 

ui>0,..., u„>0 
Ul-\ \-u n <l 

(5.1) 

n 

A random vector (Xi, ...,X„) G R" such that X, > Vi and JT, Xj =1 is 

i=l 

said to be Dirichlet distributed with parameters x n ; [107] if the 

probability density function for (Xi, X„) is given by 

i / \ r(xi + -- -+ X n U-l) x-t—l r — 1 /-i \ ^ \ r ,,— 1 

x *- ( "' "" : = r( Xl )...r( Xn+1 ) ^ ■■■< (i-E«*)^ 1 - 



i=l 



= T(xi + +fn+i X! l ... u *n l u x ^ \ provided that 

u n+1 = l-ui u n . (5.2) 

To get some feeling of such defined distribution, we notice the following pecu- 
liar aspects of this distribution. For any discrete distribution, we know that 
the probability Pi must be normalized, that is Y]; Pj = 1. Thus, the Dirichlet 
distribution is dealing with averaging of the probabilities! Or, better, it is deal- 
ing with the problem of effectively selecting the most optimal probability. The 
most primitive of these probabilities is the binomial probability given by 

p m =(^)p m (l-p) n - m , m = 0,1, 2, ....,». (5.3) 

If X is the random variable subject to this law of probability then, the expec- 
tation E(X) is calculated as 

n 

E(X) = 2J mp m = np = /i. (5.4) 

m— 1 



l As discussed recently in our work [25] on the Poincare' and geometrization conjectures. 
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Consider such a distribution in the limit: n — » oo. In this limit, if we write 
p = fi/n , then the Poisson distribution is obtained as 

Vra = ^e^. (5.5) 
to! 

Next, we notice that ml = T(m + 1). Furthermore, we replace to by real valued 
variable a and fi by x. This allows us to introduce the gamma distribution with 
exponent a whose probability density is 

Px(x) = -^-/-'e- (5.6) 

for some gamma distributed random variable X. Based on these results, we 
would like to demonstrate now how the Dirichlet distribution can be represented 
through gamma distributions. Since the gamma distribution originates from the 
Poisson distribution, sometimes in literature the Dirichlet distribution is called 
the Poisson-Dirichlet (P-D) distribution [108]. To demonstrate the connection 
between the Dirichlet and gamma distributions is relatively easy. Following 
Kingman [108], consider a set of positive independent gamma distributed ran- 
dom variables Y\, ...,Y n+ i with exponents oc\, Furthermore, consider 
the sum Y = Yi + ■ ■ ■ + Y n+ i and construct a vector u with components: 
Ui = y-. Then, since X^TJl 1 u i = 1> the components of this vector are Dirichlet 
distributed and, in fact, independent of Y. Details are given in Appendix A. 

Such described Dirichlet distribution is an equilibrium measure in various 
fields ranging from spin glasses to computer science, from linguistics to genetics, 
from forensic science to economics, etc. [103-105]. Furtheremore, most of 
fragmentation and coagulation processes involve the P-D distribution as their 
equilibrium measure. Some applications of general theory of these processes 
to nuclear and particle physics were initiated in already mentioned series of 
papers by Mekjian, e.g. see [106]. To avoid duplications, we would like to 
rederive some particular results of Mekjian's differently in order to exibit their 
connections with the previous section. 

5.3 The Ewens sampling formula and Veneziano ampli- 
tudes 

This formula was discussed by Mekjian in [109] without any reference to the 
P-D distribution. It is discussed in many other places, including Ewens own 
monograph [110]. Our exposition follows work by Watterson [111] where he 
considers a simple P-D average of monomials of the type generated by the 
individual terms in the expansion^! 

u n = ( Ul + --- + u k ) n = V ft*? 1 (5-7) 

n={n 1 ,...,n k ) 

22 Very recently Watterson's results were successfully applied to some problems in economics 
[112]. 
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This type of expansion was used in our work [83] (equations (2. 9), (2. 11)) for 
calculation of multiparticle Vencziano amplitudes. Not surprisingly, Watter- 
son's calculation also results in the multiparticle Veneziano amplitude. Upon 
multiplication by some combinatorial factor in a well defined limit such an am- 
plitude produces the Ewens sampling formula playing a major role in genetics. 
Although in Appendix B we reproduce the Ewens sampling formula (equation 
(E.6)) without use of the P-D distribution, Kingman [113] demonstrated that 
"A sequence of populations has the Evens sampling property if and only if it 
has the P-D limit" That is to say, the Ewens sampling formula implies the P-D 
distribution and vice versa. In the context of high energy physics it is the same 
as to say that the law of conservation of energy-momentum which must hold 
for any scattering amplitude leads to the P-D distribution or, equivalently, to 
the Veneziano-type formula for multiparticle amplitudes. Hence, we expect that 
our readers will consult Appendix B prior to reading of what follows. Further- 
more, since the vector u is P-D distributed, it is appropriate to mention at this 
point that equation (5.7) represents genetically the Hardy- Weinberg law [110] 
for mating specieJ^l. Hence, the Ewens sampling formula provides a refinement 
of this law accounting for mutations. 

Considear a special case of (5.2) for which X\ = = ■ ■ -xk+i — £ and let 
e = 9/ K with parameter 9 to be defined later. Then, (5.2) is converted to 

T((K -)- 1)£") 

Pxx,...,x. K ( u i, —>uk) = fc (u) = K+1 Y[ u 6 - 1 , provided that 

\ £ )\ i=l 

1 = V ^ (5.8) 
— ^i—i 

In view of (5.7), let us consider an average P(n\, tik ) over the simplex A 
(defined by YliJi u i = 1) given by 

P(ni, ...,njc) = — ; — r 1 r / " ' ( U T ' ' ' u K K 4>k( u )dui ■ ■ ■ du K - (5.9) 

ni\n 2 \ ■ ■ ■ n K \ J J 
A 

A straightforward calculation produces: 

rfr>~,n*) = nM ... nKl [T{£)f+1 U. r{{K + 1)e + n y ( 5 - 10 ) 

Up to a prefactor, the obtained product coincides with the multiparticle Veneziano 
amplitude discussed in our work [83]. To obtain the Ewens sampling formula 
(equation (B.6)) from (5.10) a few additional steps are required. These are: a) 
we have to let K — + oo while allowing many of n^s in (5.7) to become zero (this 
explains the meaning of the word "sampling"), b) we have to order remaining 
n^s in such a way that nn\ > n^ 2 ) ^ • ' ' ^ n (k) > 0,0, ...,0, c) we have to 
cyclically order the remaining n^s in a way explained in Appendix B by intro- 
ducing c'j^s as numbers of remaining n'^s which are equal to i. That is we have 



'E.g. see Wikipedia where it is known as Hardy- Weinberg principle. 
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to make a choice between representing r = Yli=i n (i) or r = Y^i=i under 
condition that k = 2l=i c i-> d) finally, just like in the case of Bose (Fermi) 
statistics, we have to multipy the r.h.s.of (5.10) by the obviously looking combi- 
natorial factor M — K\/[(c\\- ■ -c r l)((K — fc)!]. Under such conditions, we obtain: 

T({K + l)e) ~ I\0), T({K + l)e + r) = T{6 + r), r(e +" w) = Less trivial 

is the result: K\/\(K - k)\ [T(s)f] -> 6 k . Evidently, the factor — — — 

n\\ni\ ■ ■ ■ tik'- 

rl 

in (5.10) now should be replaced by . Finally, a moment of thought 

r-l n(1) ' ' ' 

causes us to replace n'^s by i Ci |ff| in order to arrive at the Ewens sampling 
formula: 

rl r 9 Ci 

P(k-n {1) ,...,n ik) ) = w Y[— ? (5.11) 
1 ' i=i % ' 

in agreement with (B.6). This derivation was made without any reference to 
genetics and is completely model- independent. To demonstrate connections 
with high energy physics in general and with Veneziano amplitudes in particular, 
we would like to explain the rationale behind this formula using an absolute 
minimum of facts from genetics. 

Genetic information is stored in genes. These are some segments (locuses) 
of the double stranded DNA molecule. This fact allows us to think about the 
DNA molecule as a world line for mesons made of a pair of quarks. Phe- 
nomenologically, the DNA is essentially the chromosome. Humans and many 
other species are diploids. This means that they need for their reproduction 
(meiosis) two sets of chromosomes-one from each parent. Hence, we can think 
of meiosis as a process analogous to the meson-meson scattering. We would 
like to depict this process graphically to emphasize the analogy. Before doing 
so we need to make a few remarks. First, the life cycle for diploids is rather 
bizarre. Each cell of a grown organism contains 2 sets of chromosomes. Maiting, 
however, requires this rule to be changed. The gametes (sex cells) from each 
parent carry only one set of chromosomes (that is, such cells are haploid !). The 
existence of 2 sets of chromosomes makes individual organism unique because 
of the following. Consider, for instance, a specific trait, e.g. "tall" vs "short". 
Genetically this property in encoded in some gen^fl- A particular realization 
of the gene (causing the organism to be, say, tall) is called "allele" . Typically, 
there are 2 alleles -one for each of the chromosomes in the two chromosome set. 
For instance, T and t (for "tall" and "short"), or T and T or t and t or, finally, 
t and T (sometimes order matters) . Then, if father donates 50% of T cells and 
50% of t cells and mother does same, the offspring is likely going to have either 
TT composition with probability 1/4, or tt (with probability 1/4) or tT (with 
probability 1/4) and, finally, tt with probability 1/4. But, one of the alleles is 
usually dominant (say, T) so that we will see 3/4 of tall people in the offspring 
and 1/4 short. What we just described is the essence of the Hardy- Weinberg 

24 This is so because the Ci numbers count how many of n'^s are equal to i. 
25 Or in many genes, but we talk about a given gene for the sake of argument. 
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law based, of course, on the original works by Mendel. Details can be found in 
genetics literature [110]. 

Let us concentrate our attention on a particular locus so that the genetic 
character(trait) of a particular individual is described by specifying its two genes 
at that locus. For TV individuals in the population there are 2N chromosomes 
containing such a locus. For each allele, one is interested in knowing the propor- 
tion of 2N chromosiomes at which the gene is realized as this allele. This gives a 
probability distribution over the set of possible alleles which describes a genetic 
make-up of the population (as far as we are only looking at some specific locus). 
The problem now is to model the dynamical process by which this distribution 
changes in time from generation to generation accounting for mutations and 
selection (caused by the environment). Mutation can be caused just by change 
of one nucleotide along the DNA strancp. Normally, the mutant allele is inde- 
pendent of its parent since, once the mutation takes place, it is very unlikely 
that the corrupt message means anything at all. Hence, the mutant can be 
either "good" (fit) or "bad" (unfit) for life and its contribution can be ignored. 
If u is the probability of mutation per gene per generation then, the parameter 
8 = 4Nu in (6.11). With this information, we are ready to restore the rest of the 
genetic content of Watterson's paper [111]. In particular, random P-D variables 
X 1; X2, ...,~K-k denote the allele relative frequences in a population consisting 
of K alleles. Evidently, by construction, they are Dirichlet-distributed. Let 
K — > oo and let k be an experimental sample of representative frequencies 
k <C K. The composition of such a sample will be random both, because of the 
nature of the sampling process, and because the population itself is subject to 
random fluctuations. For this reason we averaged the Hardy- Weinberg distribu- 
tion (5.7) over the P-D distribution in order to arrive at the final result (6.11). 
This result is an equilibrium result. Its experimental verification can be found 
in [110]. It is of interest to arrive at it dynamically. This is accomplished in the 
next subsection but in a different context. Based on the facts just discussed, 
it should be clear that both genetics and physics of meson scattering (for which 
Veneziano had proposed his amplitude) have the same combinatorial origin. All 

random processes involving decompositions r = X)j=i n (i) ( or r = * c «) are 

the P-D processes [103-105]. 

To conclude this subsection, we would like to illustrate graphically why ge- 
netics and physics of meson scattering have many things in common. This is 
done with help of the following 3 figures. 

5.4 Stochastic models for second order chemical reaction 
kinetics involving Veneziano-like amplitudes 

The role of stochastic processes in chemical kinetics was recognized long ago. 
A nice summary is contained in the paper by McQuarrie [118]. The purpose 

26 The so called "Single Nucleotide Polymorphism" (SNP) which is detectable either elec- 
trophorctically or by DNA melting experiments, etc. 
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Figure 1: The simplest duality diagram describing meson-meson scattering [114]. 
The same picture describes "collision" of two parental DNA's during meiosis 
and can be seen directly under the electron microscope. E.g. see Fig. 2. 3 in [115], 
page 18. 



— > * * 



Figure 2: Non-planar loop Pomeron diagram for meson-meson scattering [116]. 
The same diagram describes homologous DNA recombination, e.g. see Fig. 2. 2 
in [115], page 17. 
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Figure 3: The planar loop meson-baryon scattering duality diagram. The same 
diagram describes the interaction (scattering) between the triple and double 
stranded DNA helices [117]. 

of this subsection is to connect the results in chemical kinetics with those in 
genetics in order to reproduce the Veneziano (or Veneziano-like) amplitudes as 
an equilibrium measures for the underlying chemical/biological processes. Fol- 

lowing Darvey et al [119] we consider a chemical reaction A + B ^ C + D 

k-i 

analogous to the meson-meson scattering processes which triggered the discov- 
ery of the Veneziano amplitudes [13]. Let the respective concentrations of the 
reagents be a, b, c and d. Then, according to rules of chemical kinetics, we obtain 
the following "equation of motion" 

^ = -k 1 ab + fc_ic<£ (5.12) 

This equation has to be supplemented with the initial condition. It is obtained 
by accounting for mass conservation. Specifically, let the initial concentrations 
of reagents be respectively: a = A(0),/3 = B(0),j = C(0) and 5 = D(0). 
Then, evidently, a + f3 + 7 + 5 = a + b + c + d, provided that for all times 
a>0,6>0, c>0 and d > 0. Accounting for these facts, equation (5.12) can 
be rewritten as 

^ = (fc_i-fci)a 2 -[fci(/3-a) + fc_ 1 (2a + 7 + (5)]a + fc_ 1 (a + 7 )(a + /3). (5.13) 

Thus far, this is a standard result of chemical kinetics. The new element emerges 
when one claims that the variables a, b, c and d are random but are still subject 
to mass conservation. Then, as we know already from previous subsections, 
we are dealing with the P-D-type process. The new element now lies in the 
acknowledging the fact that this process is dynamical. Following Kingman [120] 
we would like to formulate it in precise mathematical terms. For this purpose, 
we introduce the vector p(t)=(pi(t),..., Pfc(t)) such that it moves randomly on 
the simplex A defined by 

A = {p(t); Pj > 0,5^^ = 1} (5.14) 
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In the present case the possible states of the system at time t which could lead 
to a new state specified by a,b,c,d at time t + At involving not more than one 
transformation in the time interval At are [119] 

o+l 6+1 c- 1 d- 1 

a- 1 6-1 c+1 d+1 | . (5.15) 
a 6 c d 

In writing this matrix, following [119], we have assumed that random variables 
a, 6, c and d are integers. Using (5.15) we obtain the following equation of motion 

P(a, 6, c, d; t + At) - P(a,b,c,d;t) = [h(a + l)(b + l)P(a + 1, b + 1, c - l,d- l;f) 

+fc_i(c+ l)(d+l)P(o- 1,6- l,c+l,d+ l;i) 
-(fciofe + fc_icd)P(a, 6, c, d; t)]At + G>(At 2 ()5.16) 

In view of the fact that the motion is taking place on the simplex A, it is 
sufficient to look at the stochastic dynamics of just one variable, say, a (very 
much like in the deterministic equation (5.13)). This replaces (5.16) by the 
following result: 

j t P a (t) = fci[(a + l)(o+l + /9-a)P 0+ i(t) + fc_i[(7 + a-o+l)(5 + a-o+l)P„_i(t) 

— [fcia(/3 — a + a) + fc_i(7 + a — a)(5 + a — a)]P a (t); provided that 
P«(0) = 1, a = a and P a (0) = if a ^ a. (5.17) 

To solve this equation we introduce the generating function G(x, t) via 

G(x,t) -^P a (t)x Q 

and use this function in (5.17) to obtain the following Fokker-Plank-type equa- 
tion 

^G{x,t) = xil-x^h-xk-^-^G+il-x^hip-a + l) 

+k_ 1 (2a + 1 + 5-l)x]^-G 
ox 

-k- 1 {a + i){a + 5){l-x)G{x,t). (5.18) 

This equation admits separation of variables: G(x, t) = S(x)T(t) with solution 
for T(t) in the expected form: T(t) = exp(—X n kit), leading to the equation for 
S(x) 

x(l-x)(l-Kx)^S(x)^P-a+l+K(2a+^+S-l)x](l-x)^S-[K(a+-f){a+p){l-s)-X n ]S(x) = 

(5.19) 

This equation is of Lame-type as discussed in [44,47] and, therefore, its solution 
should be a polynomial in x of degree at most zu, where zu should be equal to 
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the minimum of (a + 7, a + 8, (3 + 7, 8 + 8). As in quantum mechanics, this 
implies that the spectrum of eigenvalues A„ is discrete, finite and, a priory 
nondegenerate. Among these eigenvalues there must be Ao = since such an 
eigenvalue corresponds to the time-independent solution of (5.19) typical for the 
true equilibrium. Hence, for this case we obtain, instead of (5.19), the following 
final result: 

d 2 d 
x{l-Kx)-^S{x) + [f3-a+l+K{2a+ 1 +8-l)x] — S-[K{a+ 1 ){a+l3)]S = 0, 

(5.20) 

where K = k-i/k\. This constant can be eliminated from (5.20) if we rescale 
x : x — ► Kx. After this, equation (5.20) acquires the already familiar (e.g. see 
(4.41)) hypergeometric form 

d 2 d 
x(l-x)—^S(x) + U3-a+l + (2a+j+8-l)x]— S(x)-(a+~/)(a+p)S(x) = 0. 
ax z dx 

(5.21) 

In [120] Kingman obtained the Fokkcr-Planck type equation analogous to our 
(5.18) describing the dynamical process whose stable equilibrium is described by 
(5.21) (naturally, with different coefficients) and leads to the P-D distribution 
(5.2) essential for obtaining the Ewens sampling formula. Instead of reproducing 
his results in this work, we would like to connect them with results of our Section 
4. For this purpose, we begin with the following observation. 

5.4.1 Quantum mechanics, hypergeometric functions and P-D dis- 
tribution 

In our works [37,46] we provided detailed explanation of the fact that all exactly 
solvable 2-body quantum mechanical problems involve different kinds of special 
functions obtainable from the Gauss hypergeometric funcftion whose integral 
representation is given by 

1 

F(a,b,c;z)= J ^(l - ^"^(l - zt)~ a dt. (5.22) 



As is well known from quantum mechanics, in the case of a disctete spectrum 
all quantum mechanical problems involve orthogonal polynomials. The ques- 
tion then arises: under what conditions on coefficients (a, b and c) can the 
infinite hypergeometric series whose integral representation is given by (5.22) 
be reduced to a finite polynomial? This happens, for instance, if we im- 
pose the quantization condition: —a = 0, 1,2, .... In such a case we can write 
(1 - zt)- a = J27=i(T a )(- l Y( zt y and use this finite expansion in (5.22). In 
view of (5.2) we obtain the convergent generating function for the Dirichlet 
distribution (5.2). Hence, all known quantum mechanical problems involving 
discrete spectrum are examples of the P-D stochasic processei^l For hyper- 
geometric functions of multiple arguments this was demonstrated in Section 

27 Evidently, in the case of continuum spectrum we are also dealing with the P-D processes 
but the corresponding hypergeometric series conain now countable infinity of terms. 
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4. Thus, all quantum mechanical, quantum field-theoretic and string-theoretic 
processes are the P-D processes. As such they fall into a much larger class of 
stochastic processes known as random coagulation and fragmentation 
processes. We would like to conclude this section with the following additional 
observations. 



5.4.2 Hypergeometric functions, Kummer series expansions and Veneziano 
amplitudes 

In view of just introduced quantization condition, the question arises: is this the 
only condition reducing the hypergeometric function to a polynomial ? More 
broadly: what conditions on coefficients a, b and c should be imposed so that 
the function F(a, b, c; z) becomes a polynomial's The answer to this question 
was provided by Kummer in the first half of 19th century [101]. We would 
like to summarize his results and to connect them with results presented above. 
Additional details can be found in our paper [47] . By doing so we shall reobtain 
Veneziano amplitudes for chemical process described by equation (5.21). 

According to general theory of hypergeometric equations [101], the infinite 
series for a hypergeometric function degenerates to a polynomial if one of the 
numbers 

a, b, c — a or c — b (5.23) 

is an integer. This condition is equivalent to the condition that, at least one 
of the eight numbers ±(c — 1) ± (a — b) ± (a + b — c) is an odd number. This 
observation produces 24 solutions for the Gauss hypergeometric function (5.21) 
found by Kummer. Among these he singled out 6 (generating all 24) and among 
these 6 he established that every 3 of them are related to each other. 

Let us denote these 6 functions (solutions) as u\,...,ue . Then, we can 
represent, say, ui and uq using u\ and u§ as the basis set. We can do the same 
with u\ and u§ by representing them through u-i and ug and, finally, we can 
connect uj, and with u\ and u§. Hence, for our purposes, it is sufficient to 
consider, say, ui and u§. We obtain, 



u e ) \ M\ Ml 



(5.24) 



with M\ = na + b-c+l)T(l-c) = r(a + 6 + l- C )r(c-l) = 

1 T{a + l-c)^(b-c + ^) , 2 r( a )r(6) ' 1 

T(c+l-a-b)T(l-c) 9 T(c+l-a-b)T(c-l) , , 

— — : , V r -:M% = N / \ '-. The determinant of 

r(l-o)r(l-6) ' 2 T(c-a)T(c-b) 

this matrix becomes zero if either two rows or two columns become the samd 29 !. 



28 Incidentally, in the case of K-Z-type hypergeometric equations such a problem was solved 
only in 2007 [121]! 

29 The condition for the determinant to become zero is the resonance condition. It is of 
central importance in both quntum mechanics and string theory [47]. 
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For instance, we obtain: 



r(a)r(6) _ r(a - c + i)r(& - c + 1) r(c - o)r(c - ft) _ r(i - o)r(i - b) 
r(c-i) ~ r(i - c ) an r(c-i) ~ r(i - c ) ' 

(5.25) 

For c = 1 we obtain an identity. From [119] we find that (5.21) admits 2 
independent solutions: 

. . . either F(— a — 7, — a — S, j3 — a + l;Kx), for j3 > a . . 

S(x) = { —5 . (5.26) 

or (Kx) a tJ F{-j3-"f,-/3-5,a-f3 + l;Kx), for [3 < a 

Hence, the condition c = 1 in (5.25) causes two solutions for S(x) to degenerate 
into one polynomial solution, provided that we make an identification: (3 — a 
in (5.26). Notice that to obtain this result there is no need to impose an extra 
condition: a — ?o (which, in our case, is the same as 7 = 5). 

This makes sence physically both in chemistry and in high energy physics. 
In the case of high energy physics, if the Veneziano amplitudes are used for 
description of, say, nn scattering, in [83], page 54, it is demonstrated that pro- 
cesses for which "concentrations "a = b cause this amplitude to vanish. The 
Veneziano condition: a + b + c = — 1( equation (1.5) of [83]) has its analog in 
chemistry where it plays the same role, e.g. of mass conservation. In the present 
case we have a + [3 + j + 5 = const, and the Veneziano- like amplitude obtainable 
from (5. 25), (5. 26) is given now by 

KM) = r( " Q - f { - a ~ S) U- (5-27) 
-cT(-c) 

In view of known symmetry of the hypergeometric function: F(a, b, c; x) — 
F(b,a,c;x), we also have: V c (b,a) — V c (a,b). This is compatible with the 
symmetry for Veneziano amplitude. To make the analogy with Veneziano 
amplitudes complete, we have to select the following options in (5.27): a) 
a = 0, 7 = 1,6 = 1,2,. ..or; b) a = 1, 7 = 0, 6 = 0,1,2,... These conditions 
are exactly compatible with those in (1.19) of [83] for Veneziano amplitudes. 
Finally, in view of (5.22), these are the quantization conditions needed for res- 
onances to exist as required. 



6 Conclusions 

In Einstein's time to create a unified field theory was the ultimate goal of physics. 
It should be noted though, that, for instance, Pauli was not sharing Einstein's 
beleif in the unified field theory. Abdus Salam in his Nobel Lecture given in 
1979^3 on pages 518-519 writes: " I must admit I was taken aback by Pauli's 
fierce prejudice against universalism-against what we would today call unifica- 
tion of basic forces-but I did not take it too seriously. I felt this was a legacy of 

30 Here a and b have the same meaning as in (5.22) and should not be confused with con- 
centrations. 

31 Readily available at the Google search database. 
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the exasperation which Pauli had always felt at Einstein's somewhat formalis- 
tic attempts at unifying gravity with electromagnetism- forces which in Pauli's 
phrase "cannot be joined-for God hath rent them asunder". Based on the re- 
sults discussed in this paper, we sincerely hope, that should Pauli be alive, 
and be aware of the enormous universality of the stochastic theory of random 
fragmentation-coagulation processes, he would change his position on univer- 
salism. By writing an obituary for Pauli, Landau, in fact, had inagurated this 
theory 

The idea of a universal theory does not exclusively belong to physicists 
though. We mentioned already in Section 4 that Landau's paper [67] made 
a big impact on the theory of singularities and catastrophes. Rene Thom-the 
founding father of the theory of catastrophes-also was looking for the theory of 
everything. In his book [122] he made an attempt to apply the theory of sin- 
gularities to biology. Later on Robert Gilmore using Thorn's results developed 
many applications of catastrophe theory to other disciplines [123]. It is fair to 
say that the mathematical theory of catastrophes is related to physical theory 
of random fragmentation-coagulation processes as thermodynamics is related to 
statistical mechanics. We believe, that it is this stochastic theory which really 
deserves to be called the universal theory of everything. Hopefully, without 
causing a global catastrophe, the millennium into which we are entering is go- 
ing to be the millennium of incredibly versatile succession of hallmarks in the 
theory of universal catastrophes. 

Appendix A. Connections between the gamma and Dirichlet dis- 
tributions 



Using results of our paper [83] especially, equation (3.27) of this reference, 
such a connection can be easily established. Indeed, consider n+1 independently 
distributed random gamma variables with exponents a\, a n+1 . The joint 
probability density for such variables is given by 

PYi. -V„+i(*l. ...,* n +l) = ■ ■ ■ -FT 1 V S 1 1_1 ' ' ' *n+l 1_1 - ( AA ) 

T(ai) r(a„+i) 

Let now Sj = tit, where tj are chosen in such a way that Ylt\U = 1- Then, 
using such a substitution in (A.l) we obtain at once: 

oo 

p Ul ,... Un+1 (t 1 ,...,t n+1 ) = [[ r-V*]— !— ■ ■ ■ 1 if 1 ' 1 ■ ■ ■ Ck" 1 , provided that 

J l(«i) I(a n+ i) 



1 



o 

,n+l 



Since a = a\ H — • + a n+ i, we obtain: J t a ~ 1 e~ t = T(a\ + • • • + a n+ i), so that 

o 

the density of probability (A. 2) is indeed of Dirichlet-type given by (5.2). 
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Appendix B. Some facts from combinatorics of the symmetric 
group S n 

Suppose we have a finite set X. Vx G X consider a bijection X — >X made of 
some permutation sequence: x, ir(x),ir 2 (x), ...Because the set is finite, we must 
have n m (x) = x for some m > 1. A sequence (x, tt(x),it 2 (x), .., 7r m_1 (x)) = C m 
is called a cycle of length m. The set X can be subdivided into a disjoint 
product of cycles so that any permutation tt is just a product of these cycles. 
Normally such a product is not uniquely defined. To make it uniquely defined, 
we have to assume that the set X is ordered according to a certain rule. The, 
standard cycle representation can be constructed by requiring that: a) each cycle 
is written with its largest element first, and b) the cycles are written in increasing 
order of their respective largest elements. Let N be some integer and consider 
a decomposition of N as N — ^2^ =1 rii . We say that n =(no, nx) is the 
partition of N (or n h N). The same result can be achieved if, instead we would 
consider the following decomposition of N: N — ^2=1 where, according to 
our conventions, we have Cj = c^tt) as the number of cycles of length i. The 
total number of cycles then is given by K=^ i=1 a . Define a number S(N, K) 
as the number of permutations of X with exactly K cycles. Then, the Stirling 
number of the first kind can be defined as S{N,K) := (-1) N ~ K S(N, K). The 
numbers S(N,K) can be obtained recursively using the following recurrence 
relation: 

S(N,K) = (N-1)S(N-1,K)+S(N-1,K-1), N,K>1; 5(0,0) = 1. (B.l) 
Use of this recurrence allows us to obtain the following important result: 

N 

S{N, K)x K = x(x + l)(x + 2) • • • (x + N — 1). (B.2) 

K=0 

Let now x = 1 in (B.2), then we can define the probability p(K; N) = S(N,K)/N\ 
Furthermore, one can define yet another probability by introducing a notation 
[a;]^ = x(x + l)(x + 2) ■ ■ ■ (x + N - 1). Then, we obtain: 

JV K N 

k=o m K=0 

Such defined probability Pk{N; x) can be further rewritten in view of the famous 
result by Cauchy. To obtain his result, we introduce the generating function 

jr£(x) = y ^!__^ x f . . . x c « (B.4a) 
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and require that S(N, K)x = JFj^(x). This can happen only if 

N 

L c iCi!2 c 2C2! ■ ■ ■ N c *c N \ 
Thus, we obtain: 

S(N, K) = TT t^-t , provided that X = c, and Y = ic. (B.5) 

- LJ - l Ci Ci\ ^— 'i=l ^— 'i=l 

In these notations the Ewens sampling formula acquires the following canonical 
form: 

x k JL TV! 

Pk(N; x) = , provided that K = ) Cj and N = } ic. 

\x\ i Ci cA ^— 'j=i ^— 'i=l 

L 1 i—l 

(B.6) 
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